Quantum size effects in adatom island decay 
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The decay of hexagonal Ag adatom islands on top of larger Ag adatom islands on a Ag(lll) 
surface is followed by a fast scanning tunneling microscope. Islands do not always show the expected 
increase in decay rate with decreasing island size. Rather, distinct quantum size effects are observed 
where the decay rate decreases significantly for islands with diameters of 6 nm, 9.3 nm, 12.6 nm, 
and 15.6 nm. We show that electron confinement of the surface state electrons is responsible for 
this enhancement of the detachment barrier for adatoms from the island edge. 
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Thin metal films and small quantum dots have a wide 
range of applications in modern technologies. Quantum 
size effects are found to occur when the size of the nano- 
structure becomes comparable to the de Broglie wave- 
length of electrons confined within it. Such confinement 
effects have profound implications on various nanoscale 
physical properties. 

The first observation of quantum size effects in indi- 
vidual islands was made for Pb islands on Cu(lll) 0, 
where it was shown that preferred heights occur for is- 
lands without quantum well states close to the Fermi level 
0. Also Ag and Pb thin films on semiconductor sub- 
strates showed magic heights 0, 0, @ . The detection 
of quantum states was extended towards thinner films 
and a correlation between the electronic property and the 
thickness of an individual island was demonstrated . It 
has been proposed that the competition between quan- 
tum confinement, charge spilling and interface-induced 
Friedel oscillation defines the existence of the character- 
istic magic thicknesses ||. Thus, in a so called "electronic 
growth" mode for metals on semiconductors the energy 
contribution of quantized electrons confined in the metal 
overlayer determines the morphology of the growing film, 
and this energy contribution prevails over the strain en- 
ergy. This insight was used to tune the surface reactivity 
of magnesium films towards oxidation 0]. It was also 
observed that annealing of the Pb overlayer on Cu(lll) 
flf)j | and the Ag overlayer on Fe(100) ^3 leads to magic 
layer thicknesses associated with quantum size effects. 
All these observations originate from the confinement of 
electrons normal to the surface and it was implicitly as- 
sumed that the extension of the macrostructure was too 
large in the lateral direction to produce lateral confine- 
ment. However, electron confinement of the surface state 
electrons in the lateral dimension leads to quantum in- 
terference patterns on nanoscale islands. These patterns 
were observed in both homoepitaxial |l^ and heteroepi- 
taxial systems. 

During the initial growth of adatom islands the system 
is out of equilibrium since islands nucleate and since it 
cost free energy to form the step edges. When nucleation 



and growth are terminated, atoms flow from steps with 
high curvature to steps with low curvature making small 
islands shrink, while large islands grow at their expense. 
This coarsening process, referred to as Ostwald ripening, 
is a general feature and also occurs between islands in 
different layers. The adatom motion between islands of 
different layers is different from adatom motion for is- 
lands within the same layer due to the so-called Enrich- 
Schwoebel (ES) barrier |b|. Atoms that detach from the 
top island have to overcome this additional energy bar- 
rier before they can be incorporated into the bottom is- 
land. While the understanding of Ostwald-ripening pro- 
cesses between islands in the same layer has advanced 
significantly in recent years for metals [TEl as well 
as for semiconductors |Tj|, the ripening between islands 
on different layers has not been investigated thoroughly, 
although post-annealing of three-dimensional films is of- 
ten used as a simple way of forming smooth surface lay- 
ers. Only one specific aspect of the ripening between 
layers has attracted a lot of interest since it has 

been shown that on Cu(lll) and Ag(lll) surfaces, ripen- 
ing increases by orders of magnitude when the distance 
between the island edg es falls below a critical width 
El El El EJ On Cu(lll) it is found that 
this critical distance coincides with the distance where 
the surface state is depopulated. For Ag(lll), however, 
the critical distance for the occupation of the surface 
state is far above the experimentally determined criti- 
cal width Thus, these observations cannot 
be explained by quantum size effect and the underlying 
mechanism is still under intense debate |2^j . 

In this Letter we show how lateral electron confinement 
influences the decay of individual adatom islands, result- 
ing in distinct quantum size effects. We have studied in 
detail the decay of islands adsorbed concentrically on top 
of larger islands. By means of fast-scanning STM, we find 
a transition between attachment-limited decay for large 
islands, where the decay is determined by the interface to 
the bottom island, to diffusion-limited decay for smaller 
islands, where the decay is dominated by the ability of 
the adatom to diffuse to the next sink. The shift between 
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these two limiting regimes can be well understood within 
an ordinary continuum model and is caused by the in- 
creasing distance between the two step edges of the top 
and bottom islands during the decay. Very surprisingly, 
we find that the island decay is reduced significantly at 
certain island sizes. We relate these distinct island sizes 
to the absence of quantum well states near the Fermi level 
and conclude that the energy barrier for adatom detach- 
ment depends on the quantum confinement of electrons. 
Thus, we present the first observation of the influence of 
electron confinement onto surface kinetics. 

The experiments were performed on Ag(lll) single 
crystal surfaces, which were cleaned by IkeV Ar + sput- 
tering for 30 min and annealing to ~ 650° C for 30 min. 
Several sputter-anneal cycles were followed by a final 
flash to ~ 850° C, leading to terraces with a width larger 
than 1000 nm. Up to 7 monolayers (ML) of silver was 
deposited from an evaporator from Oxford Instruments 
with a rate between 0.3 and 0.8 ML/min resulting in 
stacks of adatom islands (25|. The subsequent change 
in morphology of the islands during annealing experi- 
ments was recorded by means of the home-built variable- 
temperature Aarhus STM and visualized in the form of 
so-called STM movies, i.e. sequences of time-lapsed STM 
images [jfij . As in earlier similar studies special precau- 
tion was taken to avoid any influence of the scanning 
STM tip 0,H3|. The microscope is operated in a stan- 
dard ultra-high vacuum system (base pressures of 5TQ 
Mar) that was equipped with standard surface science 
techniques for sample preparation and characterization. 
The high thermal stability of the STM combined with a 
special, active thermal drift-compensated routine allowed 
us to follow dynamically the morphological changes on 
surfaces over extended time periods (> 13h) [l^ Eil Ef^ . 

After obtaining stable imaging conditions, the dynamic 
coarsening of the islands was followed by repeatedly scan- 
ning the same spot of the surface. Figure^, shows the 
decay of the top-layer of a four-layered Ag island at 330 
K. This particular island shows the decay characteristics 
expected for island decay in the diffusion limit: The de- 
cay rate smoothly increases with decreasing island size. 
Often, however, the decay curves exhibit regions in which 
the decay proceeds at a slower rate (Fig. ^ and c) . This 
decrease in decay rates occurs at specific island sizes, e.g. 
for the island in Fig. \I]p at » 80 nm 2 . In this example 
the lower adatom island gets smaller during the island 
decay. Figure demonstrates that the decreased decay 
rates also occur for a lower island, whose area is increas- 
ing. 

Before discussing the decreased decay rates, we will 
discuss the shape of the decay curves aside from the 
plateaus. The decay of adatom islands on terraces [TilE^ I 
follows a functional dependence of the island area A on 
time of A = a ■ (to — t) 2b . The decay exponent b, ob- 
tained by fitting this power law to the island in Fig. 
varies continuously from b — (0.27 ± 0.01) for the total 
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FIG. 1: Decay of top-layer island of an adatom island stack: 
snapshots of STM movie and area development (a) 330 K, 
initial size of bottom island 2500 nm 2 , 40 pA, 2 V (b) RT, 
initial island size of bottom island: 3500 nm 2 , 5 nA, 0.3 V (c) 
RT, island size of bottom island changes from 2800 nm 2 to 
4000 nm 2 during top-layer island decay, 1.5 nA, -IV. 



curve over (0.29 ± 0.02) from t > 1000 s to (0.3 ± 0.1) 
from t > 1500 s. The last part of the decay curve for 
the island in Fig. \T]p has a decay exponent, which in- 
creases continuously from 0.2 (considering the last 1500 
s) to ks 0.25 for the last 1000 s to > 0.3 for the last 500 
s. This variation in the decay exponent with island size 
is in variance with previously observed fixed exponents 
for island dec ay o n terraces or within vacancy islands 

m e m m s m m m . 

We account for the varying decay exponent by solv- 
ing the diffusion equation in an ordinary continuum de- 
scription, which has already been used with success to 
describe the decay of adatom and vacancy islands within 
large vacancy islands and on large terraces 0, |2^, [2^ . 



3 



(a) 



_r 



>—> 



r 

i 1 



1_ 



R 




FIG. 2: (a) Theoretically considered scenario, for details see 
text (b) Decay curve from integration for the island in Fig. 
la (dashed line) with the value of prefactors as indicated (full 
lines); integration values: Ar — 2560 nm 2 , E s = 0.13 eV, 
E e = 0.71 eV, 7 = 0.75/nm, T = 330 K. 



In this model, we consider the geometry sketched in Fig. 
|2K [25| : A small adatom island of radius r is placed con- 
centrically on top of a larger adatom island of radius R. 
The decay of the top-layer island can be described by 
the three net fluxes indicated in Fig. (i) The flux of 
adatoms from the adatom island to the terrace J r , (ii) 
the diffusion flux of adatoms over the terrace J r ^R, and 
(iii) the attachment flux of adatoms to the outer island 
Jr. Considering mass conservation and equality of net 
fluxes at steady state in analogy to the theoretical mod- 
eling of a vacancy island within a vacancy island [lf)j . 
yields a differential equation for the island area A = irr 2 : 



dA 
~dt 



2nD p eq {r) - Peq{R) 
n P r + P r ^R + Pr 



(1) 



where n is the atomic density in the surface layer, D is the 
diffusivity of adatoms over the surface, and p e q(r) = Poo • 
exp( k ^ nf ) is the equilibrium adatom concentration of a 
step of curvature 1/r with k the Boltzmann constant and 
T the absolute temperature. The prefactors P r = a/r, 
Pr-iR = ln(R/r), and Pr = a/Rs result from the three 
net fluxes involved (Fig. EH) with a the surface lattice 
constant and s — sq ■ e ~ E ^/ kT ; with E s the ES barrier. 

P r is by far the smallest term for our experimental con- 
figuration, reaching at most a tenth of Pr, and it is there- 
fore negligible (E s = 0.13 eV, s = 0.25 for Ag(lll)) 0. 
The relative importance of the other two terms changes 
during the decay. A dominance of Pr is indicative of 
an attachment- limited decay, i.e., the decay kinetics is 
dominated by the ability of atoms to attach to the bot- 
tom island. If, however, P r ^R. dominates, the decay is 
diffusion- limited, i.e., the decay kinetics is dominated by 
the diffusion of the adatoms between the island edges. 

We can linearize the exponential expressions for the 



300 
250 
200 
■f 150 
< 100 
50 



(a) 



* ** 




1000 2000 3000 4000 

150 tSL 




FIG. 3: (a) Areas of several top-layer islands decaying at 
RT (b) Histogram of measured island sizes for 10 different 
islands; island sizes have been measured every 10 s; vertical 
bars indicate island sizes, which have no quantum well state 
at the Fermi level (see text). 



adatom concentrations of eq.n an d thereby obtain an an- 
alytical solution in the attachment-limited case (Pr >> 
Pr^R.) which yields A — a ■ (to — t) 2b with an invari- 
able decay exponent of b = 1/3. However, the increasing 
importance of the diffusion-limited term for decreasing 
island size results in a mixture of the two decay regimes, 
and thus in a variety of 'apparent' decay exponents, in 
accordance with our experimental observation. A correct 
description of the experimental data in principle demands 
a numerical integration of eq.^ as is shown in Fig. [3) for 
the experiment of Fig. using the material parameters 
of Ag(lll) determined previously 0|. The increasing 
importance of the diffusion-limited term is evident from 
the increasing value of the prefactor P r ^R. 

The island size f at which the decay will be predomi- 
nantly diffusion-limited (Pr < P t ^r) is given by: 



r = R- e Rs 



(2) 



The type of island decay depends thus on the sizes of the 
two islands involved (r and R) and the Ehrlich-Schwoebel 
barrier, i.e. the material and specific crystal face. 

The continuum theory thus explains the observed con- 
tinuous variation of the decay exponent. It can not, how- 
ever, explain the observed decrease in decay rate at cer- 
tain distinct island sizes (Fig. Such apparent de- 
cay plateaus are observed for many islands for lower is- 
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land sizes varying from 150 to 10.000 nm 2 regardless of 
whether the lower island grows or decays. In Fig. [3^ 
we show some more decay curves. These reveal a dis- 
tinct reduction in the decay rate by approximately one 
order of magnitude for different islands at certain island 
sizes. Although the different islands show plateaus at 
similar values, not all decay curves show plateaus at all 
the characteristic island size values. 

In fig.|2t> we have plotted the diameter from the decay 
curves measured at fixed time intervals for several islands 
in a histogram. Distinct preferred sizes are observed for 
islands diameters of 6 nm, 9.3 nm, 12.6 nm, and 15.6 
nm. These peaks correspond to those island diameters, 
where the decay rate is much slower than expected from 
the continuum treatment. 

For the vertical growth mode reported in 2] the is- 
land heights with no quantum well states at the Fermi 
level are preferred. In the lateral dimension, magic island 
sizes without quantum well states of the surface state 
electrons near the Fermi level are given approximately 
by A F /4 + n ■ A F /2. With A F /2 = 3.7 nm for Ag(lll), 
which leads to 5.55, 9.25, 12.95, and 16.65 nm for n = 1 
... 4. The comparison to the maxima in Fig.|3jD reveals a 
clear correlation between the size of the stabilized islands 
and these values. 

Thus, we have shown that the island decay slows down, 
whenever there is no quantum well state near the Fermi 
level. This implies that not only the total energy of an 
island oscillates as a function of its size, as in " electronic 
growth" P, 0, S H S 0, 0, but also the energy 
barrier for detachment does so. 

From the comparison of the decay rates, in the 'regu- 
lar' decay regime and the decreased decay rates we can 
estimate the energy difference of this detachment bar- 
rier (under the assumption of the same prefactor) to 
(0.07 ± 0.03) eV. This is about one tenth of the total 
detachment rate, which was determined in a different de- 
cay experiment to 0.7 eV 

In conclusion, the study of the dynamics of the top- 
layer of a multilayered adatom island stack on Ag(lll) 
with a fast scanning tunneling microscope reveals two 
remarkable phenomena: A transition from attachment- 
limited to diffusion-limited decay, and a quantum size 
effect in island kinetics. Both phenomena are expected 
not to be restricted to the particular system. The former 
is important whenever the island sizes are close to fulfill- 
ing equation The latter will influence island decay for 
all surfaces with occupied surface states. 
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